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Abstract 

We formulate the Bogoliubov variational principle in a mathematical framework 
similar to the generalized Hartree-Fock theory. Then we analyze the Bogoliubov theory 
for bosonic atoms in details. We discuss heuristically why the Bogoliubov energy should 
give the first correction to the leading energy of large bosonic atoms. 
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1 Bogoliubov theory 

In this section we formulate the Bogoliubov variational principle in the same spirit of the 
generalized Hartree-Fock theory [4]. Our formulation bases on the earlier discussions in 
[EJCL9]. 
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1.1 One-body density matrices 

We start by introducing some conventional notations. Let f) be a complex separable Hilbert 
space with the inner product (., .) which is linear in the second variable and anti-linear in 
the first. Let t)N '■= <S^ m ^ be the symmetric tensor product space of iV particles and let 
J 7 = ^(f)) := ©7v =0 §n be the bosonic Fock space. 

Let £>(J-") be the space of linear bounded operators on J 7 . Any quantum mechanical state 
(state for short) p : B^F) — > C is identified with a positive semi-definite trace class operator 
P on J 7 with Tr(P) = 1 in such a way that 

p(B) = Tr (BP) for all B G ^(J 7 ). 

For example, a jmre state is a state corresponding to the one-dimensional projection 

of a unit vector f 6 J 7 , and a Gz&frs state is a state corresponding to Tr(exp(— iJ)) -1 exp(— H) 

for some Hamiltonian H : J 7 — > J 7 such that exp(— iJ) is trace class. 

The dual space t)* can be identified to f) by the anti-unitary J : f) — > f)*, 

It is convenient to introduce the generalized annihilation and creation operators on f) © f)* by 

A(f®Jg) = a(f) + a*(g), 
A*(f®Jg) = a*(f) + a(g), for all f,gei) 

where a(f) and a*(f) are the usual annihilation and creation operators. Note that if we 
denote 

then we have the conjugate relation and the canonical commutation relation (CCR) 

A*(F 1 ) = A(JF 1 ) , [A(F 1 ),A*(F 2 )\ = (F U SF 2 ) for all F 1 ,F 2 G f) © f)* 

where [X, F] = XY — YX. 

Now we can define the one-particle density matrix (1-pdm for short) r : f) © f)* — >■ f) © f)* 
of a state p by 

(Fi, TF 2 ) = p(A*(F 2 )A(F 1 )) for all Fi,F 2 G f) © fj*. 
Such a 1-pdm may be also written as 

r=f 7 a 1 (i) 

V JaJ 1 + J 7 J* y 1 ; 

where 7 :!)—>■ f) and a : f)* — )■ f) are linear bounded operators defined by 

(f,19) = P(a*(g)a(f)), (f,aJg) = p(a(g)a(f)) for all f,g G f). 

It is obvious that any 1-pdm is positive semi- definite. The following lemma expresses the 
condition T > in terms of 7 and a. Its proof is provided in the Appendix. 
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Lemma 1.1. Let F be of the form (Qp. Then T > if and only if 7 > 0, a* = JaJ and 

7>aJ(l + 7)- 1 J*a*. (2) 

Remark. The fermionic analogue of the inequality (J2J) is aa* < 7(1 — 7) jl]. We do not know 
if (j2J) can be reduced to aa* < 7(1 + 7) or not. 

Of primary physical interest are the states with finite particle number expectation. Recall 
the particle number operator 

00 

■/V := } y Nlt, N = ^2 a*(u n )a(u n ) 

N=0 n 

for any orthonormal basis {u n }^ =l for f). It is straightforward to see that if a state p has the 
1-pdm of the form (TjQ) then 

PW = Tr( 7 ). 

Hence p has finite particle number expectation if and only if 7 is trace class. 
1.2 Bogoliubov transformations 

Definition (Bogoliubov transformations). A bosonic Bogoliubov transformation is a linear 
bounded isomorphism V : f) © fy* — > \) © f)* satisfying 

JVJ = V and V*SV = S. 

These conditions ensure that the Bogoliubov transformations preserve the conjugate re- 
lation and the canonical commutation relation, namely 

A*(VFi) = AiyjFi) and [A(VFi), A*(VF 2 )] = (F h SF 2 ), WF U F 2 e fj © fj*. 

The Bogoliubov transformations form a subgroup of the isomorphisms in f) © f)*; in 
particular, if V is a Bogoliubov transformation then V -1 and V* are also Bogoliubov trans- 
formations. Note that any mapping V satisfying JVJ = V must have the form 

v ={jvj JUJ>) < 3 > 

for some linear operators U :f)— »t), V:f)*— > t). 

We say that a Bogoliubov transformation V is unitarily implementable if it is implemented 
by a unitary mapping Uy : T — > F, namely 

A{VF) = V V A(F)W V for all F e t) © f)*. (4) 

The following result determines whenever a Bogoliubov transformation is unitarily imple- 
mentable. This result is well-known and we provide its proof in the Appendix for the reader's 
convenience. For the fermionic analogue, see jl] (Theorem 2.2) . 

Theorem 1.2 (Unitarily implementable Bogoliubov transformations). A Bogoliubov trans- 
formation V:f)ffif)*— »f)ffify* of the form is unitarily implementable if and only if the 
Shale- Stinespring condition Tt^iVV*) < 00 holds. 
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Unlike to the fermionic case jl], the bosonic Bogoliubov transformations are not unitary 
mappings on f) © f)*. However, we can still use the Bogoliubov transformations to diagonalize 
some certain operators on f) © fj*. Of our particular interest is the diagonalization of the 
1-pdm's. 

Theorem 1.3 (Diagonalization 1-dpm's by Bogoliubov transformations). IfT has the form 
(CP with T > and Tr(y) < oo then for an arbitrary orthonormal basis {u n } for t), there is 
a unitarily implementable Bogolubov transformation V : f) © f)* — > f) © f)* diagonalizing Y in 
in the basis u\ © 0, u 2 © ; © Ju\, © Ju 2 , namely 



( ^ 



VTV 



\ 



A 2 



1 + Ax 



1 + A 2 



V 



(5) 



/ 



Remark. The finite-dimensional case is Theorem 9.8 in [19]. See jl] (the proof of Theorem 
2.3) for the fermionic analogue. 

To prove Theorem 11.31 we start with a simple diagonalization lemma. This is a general- 
ization to infinity dimensions of Lemma 9.6 in 



Lemma 1.4. Let A be a positive definite operator on f) ffi \f such that J A J = A and SA 
admits an eigenbasis on f) © [)*. Then for any orthonormal basis Ui,u 2 , ■■■ for f) ; there exists 
a Bogoliubov transformation V such that the operator V*AV has eigenvectors of the form 
{w n ©0}U{0© Ju n }. 

Remark. In this result the Bogoliubov transformation V needs notbe unitarily implementable. 

Proof. 1. Let {ui} be an orthonormal basis for t). We shall define the Bogoliubov transfor- 
mation V by 

V{ui © 0) = v h V(0 © JUi) = Vi, 
where {%} U {vj} is an eigenbasis of SA such that 

(i) (v h Svj) = Sij, (v^ Svj) = -6ij and (v h Svj) = for all i,j = 1, 2, . . . 

(ii) Jvj = Vj for all j = 1,2,... 

2. Let Vi be a normalized eigenvector of SA with eigenvalue Ai- Using Avi = X±Svi we 
find that 

(vi,Av!) = Xi(vi,Sv!). 

Since A is positive definite and S is Hermitian, both of Ai and (vi, Sv\) must be real and 
non-zero. Therefore, we can normalized v\ in such a way that (v\,Sv{) G {±1}- 
Defining V\ = Jv\ and using J A J = A we have that 



SAvi = SJAvi = —JSAvi = —J\\V\ = — Xiv m+1, 
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where we have used Jhat Ai is real and that JS = —SJ. Thus v% is an eigenvector of SA 
with the eigenvalue Ai = — Ai. 

Since Ai 7^ 0, Ai and Ai must be different. On the other hand, 

Ai(ui,555i) = {v x ,AX x ) = (Av x ,v x ) = X x {v x ,Sv x ). 

Thus (v x ,Sv x ) = 0. Moreover, since 

{v x ,Sv x ) = {Jv x ,SJv x ) = {Jv x ,-JSv x ) = -(Sv x ,v x ) = -(v x ,Sv x ), 

by interchanging v x and v x if necessary we can assume that (v x ,Sv x ) = 1 and (v x ,Sv x ) — — 1. 

3. Let V = Span{ui,55i} and W = {SV) 1 - = «S(F ± ). We shall show that 

fj © fj* = V © W. 

Indeed, if a <E V HW then a EV = Spem{v x ,v x } and (a, S'u) = for all v G V. Because 
(fi,iSt>i) = 1, (ui,«Sui) = —1 and (v x ,Sv x ) = 0, we must have a = 0. Thus PI W = {0}. 

On the other hand, if a G (V© W) 1 C^nlf 1 then Sa e S'(l/ ± ) n .S^) = lfn7 = 
{0}, and hence a = 0. Therefore, (F © W) 1 - = {0}. 

Moreover, since V is finite dimensional and W is closed, the direct sum space V © is 
a closed subspace of f) © fj*. Thus E) © ^* = V © W. 

4. We prove that 5^4 maps W into itself. Indeed, using V = SAV we have WA-SV = 
S(SAV) = AV. Since A is symmetric, we get AW1.V, and hence SAWASV. Thus 

C (SV) 1 = W. 

Because SA admits an eigenbasis on f) © f)* = ^©VT and SA leaves V and W invariant, 
SA also admits an eigenbasis on W. We then can restrict SA on W and conclude the desired 
result by an induction argument. □ 

Next, we show that r + |«S satisfies all assumptions on A in Lemma [1.41 

Lemma 1.5. Let T be of the form (QP with T > and Tfi^y) < 00 and let Y x := Y + |«S. 
T/ien Ti is positive definite on f) © f)*; moreover, JT X J = T x and ST X admits an eigenbasis 
on fj © [)*. 

Proof. 1. It is straightforward to check that J7TXJ7 = Tx. We now prove that Tx is positive 
definite. 

First at all, it follows from T > that 

(/ © J<7, (r + 5)/ © Jg) = (g® Jf, F(g © J f)) > 0, 
namely r + S > 0. Thus 

T x = T + ^S=^[T+(T + S)]>0. 

Next, we check that T x is injective. Assume that there exists ip G Ker(r!)\{0}. Then 
since J and F x commute, we have Jip G KerfT^yjO}. Because J leaves the subspace 
Span{(/?, Jtp) C Ker(Ti) invariant, J must have a non-trivial fixed point in this subspace. 
Thus there exists / G f)\{0} such that © Jf) = 0. Using this equation we find that 

(/© J(t/),r(/e J(tf))) = ( / , 7/ ) + t 2 ( / ,(i + 7 ) / )-t( / ,(2 7 + i)/) 

= (t-i) 2 (f,if) + (t 2 -t)\\ff<o 
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for some t < 1 and near 1 sufficiently. However, it is contrary to T > 0. Thus Ti must be 
injective. 

To see that Ti is positive definite we can introduce T\ , the unique positive semi-definite 

1/2 1/2 

square root iy on f) © f)*. Since I\ is injective, Iy is also injective, and hence 

(<^,r» = ||r} /2 w|| 2 > for all <p^0. 

2. We show that STi has an eigenbasis on f) © f)*. Although STi is not a Hermitian, we 
may associate it with the Hermitian C = F^ST^. 

We can see that C has an orthonormal eigenbasis for f) © f)*. Indeed, it is straightforward 
to see that 



C 2 



ri /2 (srs)i7 z = rj 



1/2 



.1/2 



7 
-a" 



—a 
J 7 J* 



I 



1/2 



1/2 



7 —a 
-a* J 7 J* 



if 2 + ~r\ 

1 2 



1/2 



7 

-a* 



—a 



1/2 



7 
a* 



a 
J 7 J* 



Because 7 is trace class, aa* is also trace class due to inequality ([2]). Thus {C 2 — \ V) is a 
self-adjoint Hilbert-Schmidt operator, and hence it has an orthonormal eigenbasis on f) © h*. 
Therefore, C 2 has an orthonormal eigenbasis. Note that if <p is an eigenvector of C 2 then 
dp is also an eigenvector of C 2 with the same eigenvalue. Because C maps the subspace 
Span{<£>, dp} into itself, we can diagonalize to obtain an orthonormal eigenbasis of C on this 
subspace. By induction, we get an orthonormal eigenbasis of C on I) ffi f)*. 

1 /2 

Now note that if (p is an eigenvector of C then SY^ ip is an eigenvector of STi with the 
same eigenvalue since 



sr^srfV) = srnri /z srnv? = srncv) 



l/2/ F l/2 QF l/2 



-1/2, 



1/2 | 

Moreover, because both of S and I\ are injective, ST^ maps a basis on f) © f)* to another 
basis. In particular, ST^ 2 maps an eigenbasis of C to an eigenbasis of SYi. 

Now we can prove Theorem 11.31 similarly to Theorem 2.3 in jl]). 



1/2 



□ 



Proof of Theorem \1.3[ 1. Apply Lemma rOl with A = Ti := r+|iS, we can find a Bogoliubov 
transformation V on f) © f)* such that, with respect to the orthonormal basis {u n © 0} U {0 © 

Ju n }, 



( Ai 



VTiV 



A2 + I 



A a + J 



A 2 + i 



V 



■■■ / 



which is equivalent to (jSJ). 
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We claim that in (j^D we must have A n > and J2 n < 00 • It follows from (jSJ) and 
V*TV > that A n > 0. In order to prove the boundedness J2 n A n < 00 we note that 



TS(T + S) 



7(7 + 1) — aa* 7a — aJ'-yJ* 
a*j-J-fJ*a* a*a - 77(7 + 1)J* 



is a self-adjoint trace class operator. Using the diagonal form 

VT5(r + <s)v = [vtv] s [v*(r + <s)v] 

/ A x (Ai + l) 

A 2 (Ai + l) 



\ 



-Ai(Ai + i; 



-A 2 (A X + 1) 



V 



we conclude that A n (A n + 1) < 00, which is equivalent to J2 n ^« < 00. 

2. Finally we show that the Bogoliubov transformation V constructed above is unitarily 
implementable. Assume V has the form (j3J). Then by Theorem II .2[ it suffices to prove that 
VV* is a trace class operator on f}. 

It follows from the representation flS} that the upper left block of V*TV is a positive 
semi-definite trace class operator on f). By direct computation, we can see that the upper 
left block of 



VTV 



U* J*V*J* 
V* JU*J* 



7 a 
a* I + J7J* 



U V 
JVJ JUJ* 



is 



U*jU + J*V*J*U + U*aJVJ + J*V*VJ + J*V*-fVJ. 

Because 7 is trace class, we have U*jU and J*V* r yVJ are trace class. Thus J*V*J*U 
U*aJVJ + J*V*VJ is trace class. 

Moreover, using the Cauchy-Schwarz inequality 



Tr(XY + Y*X*)\ < 2(Tr(A:A:*)) 1/2 (Tr(rr*)) 



a/2 



we find that 



00 > Tr [U*aJVJ + J*V*J*a*U + J*V*VJ] 

= Tr [(U*aJ)(VJ) + (VJ)*(U*aJ)*} + Tr(VV*) 
> -2(Tr(U*aa*U*)) 1/2 (Tr(VV*)) 1/2 + Tr(VV*). 

Note that Tr(U*aa*U*) < 00 because aa* is trace class. Thus Tr(W*) < 00. 



□ 



1.3 Quasi-free states and quadratic Hamiltonians 

Definition (Quasi-free states). A quasi-free state p is a state satisfying Wick's Theorem, 
namely 



p[i4(*i)...A(F : 



2m- 1; 



for all m > 1 



(6) 
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and 



p[A{F x )...A(F 2m )] = p[A{F^ 1) )A{F^)]...p[A{F^ 1) )A{F a( ^ ) )] (7) 



where P2 m is the set of pairings 
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{a G S 2m | <r(2j - 1) < a(2j + 1), j = 1, . . . ,m- 1, 

a(2j-l)<a(2j), j = l,...,m}. 



A crucial point is that we have one-to-one correspondence between the set of quasi-free 
states with finite particle numbers and the set of 1-pdm's. If a quasi-free state is a pure 
state, namely a one- dimensional projection on the Fock space, we call it a quasi-free pure 
state. 

Theorem 1.6 (Quasi-free states and quasi-free pure states). 

(i) Any operator V : f) © f)* — > f) © f)* of the form (Tjp satisfying T > and Tr(j) < oo is 
the 1-pdm of a quasi-free state with finite particle number expectation. 

(ii) A pure state (\&| with finite particle number expectation is a quasi-free state if and 
only if = Uy |0) for some Bogoliubov unitary mapping Uy as in (0). 

Moreover, any operator F : f) © f)* — >• fj © f)* o/ t/ie /orm (T7]j satisfying T > and 
Tr(7) < oo i/ie 1-pdm of a quasi-free pure state if and only ifTST = —T. 

Remark. The characterization of quasi-free pure states were already proved in [19] (with a 
different proof). For the fermionic analogues see [lj (Theorem 2.3 and Theorem 2.6). 

Proof, (i) Note that the set of quasi-free states is invariant under Bogoliubov unitary map- 
pings. Indeed, if the Bogoliubov transformation V is implemented by the unitary mapping 
Uy : T — > J 7 as in (J3|) and T is the 1-pdm of a quasi-free state p then V*TV is the 1-pdm of 
the quasi-free state py*rv defined by 



Therefore, due to the diagonalization result in Theorem ll.3[ it remains to show that any 
operator of the form 



where £ is a positive semi-definite trace class operator on f), is indeed the 1-pdm of some 
quasi-free state. 

Because £ is trace class, it admits an orthogonal eigenbasis {ui}^ for t) corresponding 
to eigenvalues {Aj}^ 1 . Let I = {i e N|Aj > 0}. Then we may choose G (0, oo) such that 



p v « TV (B) := p(U v 5U v ) for all B G B(F). 




(1 - exp(-ej)) 



i 



1 + X h i G I. 



(8) 



Denote ai = a(ui) for short. Let 




(9) 



S 



where n is the orthogonal projection onto the subspace Ker^^j a*Oj] . Similarly to the 
fermionic case (see Theorem 2.3 in [1]), it is straightforward to check that T is the 1-pdm 
of the state p = Tr[G] _1 G and that p is a quasi-free state. For the reader's convenience we 
provide this part of the proof in the Appendix. 

(ii) If \I> = Uy |0) for some Bogoliubov unitary mapping Uy then using U v = Uy-i and 
(HD we have AiV^F) = W V A(F)U V . Therefore, 

(V\A(F 1 )A(F 2 )...A(F n )\V) = (0|U v A(F 1 )U v UyA(F 2 )Uy...U v A(F n )Uy|0) 

= (o| ^(v- 1 ^)^^- 1 ^)...^^- 1 ^) |0) . 

It is then obvious that the state satisfies equations ([6])-([7j) in Wick's Theorem, and 

hence it is a quasi-free state. 

Reversely, suppose that the pure state (9\ is a quasi- free state with finite particle 
number expectation. Then by the first statement in Theorem II. 6[ 

(tf | = Tr[G] -1 U v GU£ 

for some Bogoliubov unitary mapping Uy and for some G given by On the other hand, 
the only rank-one operator G of the form ([9]) is the vacuum projection |0) (0| (namely £ = 0). 
Thus, up to a complex phase, $ is equal to Uy |0). 

Now we consider the 1-dpm's of quasi-free pure states. Suppose that \1/ is a quasi-free 
pure state with finite particle number expectation and its 1-dpm is T. Due to Theorem 11.61 
there is a unitarily implementable Bogoliubov transformation V such that 



v*rv 




1 



The identity TST = — T follows from 



v*i\srv = (v*tv)(v*sv)~ l (v*tv) 

\ ! / 
1 ) b \ 1 

° J J = -VTV. 

Reversely, let Y : f) © f)* t) © ff be of the form such that T > 0, Tr(7) < oo 
and TST = —T. Then by Theorem 11.61 T is the 1-dpm of a quasi-free state and there is a 
unitarily implementable Bogoliubov transformation V such that 



VTV 



£ 
1 + J£J* 



for some positive semi-definite trace class operator £ on f). The identity T^r = — r implies 
that 

'£ \ / £ \ / £ 
1 + J£J* J \ 1 + J£J* J " V 1 + J£J* 

The only solution to this equation is £ = 0. Therefore, F is the 1-dpm of a quasi-free pure 
state with finite particle number expectation. □ 
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One of the main motivation of considering the quasi-free pure states is that they minimize 
the quadratic Hamiltonians. 

Definition (Quadratic Hamiltonian) . Let A be a positive semi-definite operator on f) © f)* 
and J A J = A. The operator 

H A =Y i {F i ,AF j )A*(F i )A(F j ), 

acting on F is called a quadratic Hamiltonian corresponding to A. Here {-Fj}i>i is an 
orthonormal basis for f) © f)* (the sum is independent of the choice of {Fj}j>i). 

Remark. (i) Alternatively, we can describe H A by 

(\I/, H A ^f) = Tr[AT\s,] for all normalized vector f G J, 
where I\ is the 1-pdm of the pure state |\&) (\&|. 

(ii) The condition J A J = A is just a conventional assumption since if this condition does 
not holds then we can consider A' = \{A + J A J) which satisfies that J A' J = A' 
and, formally, 

H A , = H A +~Tr[AS\. 

This formal formula makes sense when, for example, A is trace class. 

(ii) As we shall see below, that A > is the necessary and sufficient condition such that 
H A is bounded from below. Moreover, in this case H A > 0. 

We are interested in the ground state energy of H A , 

E(H A ) := M{p(H A )\p is a state with p(N) < 00} (10) 
Theorem 1.7 (Minimizing quadratic Hamiltonians). Let A, H A and E(H A ) as above. 

(i) We have E(H A ) = mf{p(H A )\p is a quasi-free pure state}. 

(ii) If there is a unitarily implementable Bogoliubov transformation V A such that V*^AV A is 
diagonal then there is a quasi-free pure state po such that po{H A ) = E(H A ). Moreover, 
if A is positive definite then po is unique. 

(Hi) If the variational problem UP) has a minimizer then A is diagonalized by a unitar- 
ily implementable Bogoliubov transformation V A . Moreover, if V is the 1-pdm of the 
minimizer then we have 

AT = -ASI^q) [AS). 
In particular, ATS = ST A < 0. 

Remark. (i) The above statements (i) and (ii) already appeared in [19] in the finite- 
dimensional case (in this case A is always diagonalizable by Lemma II .4j) . 
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(ii) If the operator W is not self-adjoint but U 1 WU is self-adjoint for some invertible 
operator U then we can still define the projection l(_ oo )[W / ] by 

i ( -oc,o)[w] := ui^Au^wup- 1 . 

It is easy to check that the definition is independent of the choice of U. In particular, 
we can define 

l^o) [AS] := TO-Voo^AS^)- 1 ]^ 
where V^ASiVX) 1 is self-adjoint. 

Proof, (i) We show that for any state p with finite number particle expectation, there is a 
quasi- free pure state p such that p(H A ) < p(H A ). 

By Theorem II. 3| there is a unitarily implementable Bogoliubov transformation V such 
that 

r = v ( ^ iv* 
V o 1 + Jir ) 

where £ : f) — > f) is a positive semi-definite trace class operator. Thus 



p(H A ) = Ti[AT] = Tr 



V*AV 



l + J^J* 

Because V*AV commutes with J , it has the block form 



V*AV = 

where < a : f) -> f) and 6 : f)* — >• f). Thus 
P(H A ) = Tr 



a 6 
JbJ JaJ* 



a b 
JbJ JaJ* 



£ 
1 + J£J* 



= 2Tr[<] + Tr[a] > Tr[a] 
By Theorem 11.61 there is a quasi-free pure state p whose 1-pdm is 



v(jj»)v. 

It follows from the above discussion that p(H A ) < p(H A ). 

(ii) Assume that A is diagonalized by the unitarily implementable Bogoliubov transfor- 
mation V4, namely 

' d 
' A \ JdJ* 



a = v*a : T r Tst i v A 



Y d 


° 1 


A 


JdJ* ) 



where d : f) — > () is positive semi-definite. For any state p we have 

p{H A )=Tr[AT]=Tr 
where T is the 1-pdm of p. We may write V4IVV4 in the block form 
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where < 7 :()—)• f) and a : f)* — > P). Thus 
P(^) = Tr 



d 
JdJ* 



7 
a* 



a: 

1 + JlJ* 



2Tr[cfy] +Tr[d] > Tr[d] 



Denote by po the quasi-free pure state having the 1-dpm 




1 



vr 1 



sv* A 




1 



V A S. 



Then p (H A ) = Tr[d\ and hence p is a ground state of 

Moreover, if A is positive definite then Tr[c?7] > unless 7 = 0. Therefore, po is the 
unique ground state of H A among the quasi-free states. 



(iii) Assume that problem ffTOj) has a minimizer and T is the 1-dpm of the minimizer. 

1. We first prove that AS and ST commute. Let a be an arbitrary trace class operator on 
f)©f)* such that a = a* = JaJ. It is straightforward to check that exp(ieHS) is a Bogoliubov 
unitarily implementable transformation for any £6i Similarly to the variational argument 
for Hartree-Fock-Bogoliubov theory in [12] (p. 284), we consider the trial states 

T e := exp(-z£«Sa)rexp(iea«S) = T + e[iTaS - iSaT] + 0{e 2 ), eel. 

Since e = minimizes the functional e h-> Tr[*4.(T £ — T)] we find that 

= A Tr[A{T £ - T)] = Ti[aB] with B := iSAT - iTAS. 

Note that B = B* = JBJ. 

Now let b be any trace class operator on fj© f)*. Since a := 6 + 6* + JbJ + Jb* J satisfies 
that a = a* = JaJ ', we have 

= Tr[afi] = 49?Tr[65]. 

By changing 6 a complex phase, we conclude that Tr[65] = for any trace class operator b. 
This implies that B = 0. Thus AST = FAS, namely [AS, ST] = 0. 

2. Now let Vr be the unitarily implementable Bogoliubov transformation such that 

for some trace class operator £ > on t). Thus the operator V r ~ 1 5rVr = iSVpTVr leaves the 
spaces t) © and © I)* invariant . Since ASVr commutes with Vf 1 <SrVr, it also leaves 
the spaces f) © and © ()* invariant. Moreover, since V^ASVr commutes with J, it must 
have the form 

V^ASVr = ( I J dJ , ) . 
Using Vp 1 = SV^S we then conclude that 

SVlSASVrS = ( I j° dJ , ) 
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where d > since A > 0. Thus the unitarily implementable Bogoliubov transformation 
V A := SVyS diagonalizes A. 

3. Finally we prove that AT = — ASl(_ 00j o)[.AS]. Denote 



for any 1-dpm T'. We find that 



7 a 
(5)* JjJ* 



and V^r'O^* 



7 
(a'] 



a 

I T* 



< Ti[A(T' -T)] = Ti[(V* A AV A )V A \T' -T)^)- 1 } 
d 



Tr 



JdJ* 



7—7 a — a 
{a' -a)* J (7' -7)J* 



2Tr[d7 / ] -2Tr[d7]. 



Because this inequality holds true for any positive semi-definite trace class operator 7' on 
f), we conclude that Tr(cFy) = 0. By writing TrfcFy] = Tr[(d 1/2 7 1/2 )*d 1/2 7 1/2 ] we obtain 
d 1 / 2 ^/ 1 / 2 = 0, and hence oFy = 0. This also implies that da = since 



Thus 






] 




( d 


° ) 





JdJ* 




V 


JdJ* J 



(ad)* (ad) = d(a)*ad < d(l + \ \j\\c(t)))ld — 0- 

7 a 
(a)* I + J7J 

= (v a av a )(v a 1 t(v* a )- 1 ) = v* a at(v 

It can be rewritten as 

at = (V 

Moreover, since d > we find that 

d n \ 

1 





JdJ* 




JdJ* 



d 
JdJ* 



L (-oo,0) 



d 
-JdJ* 



S 



-(VyiV^lc-oco) [SV* A AV A ]S 

-V^V^l ( _oa,0) [(SV A )- 1 (AS)SV A \ S 

-V* A AV A (SV A )- 1 l ( _ oofi) (AS)SV A S 
-V^l ( _ 00 ,o)(^)(V^)- 1 . 



Thus (HU) can be rewritten as 



AT = -SAl(-oo,o)[AS]. 
Moreover, it follows from (TTTT) that 



ATS = (Vi)"' ( ° ) V^S = SV A (l J dJ , ) 



V* A S < 0. 



□ 
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1.4 Bogoliubov variational theory 

The Bogoliubov variational states should include not only the quasi-free states (like the 
Hartree-Fock theory) but also the coherent states, which correspond to the condensation. 
To describe the formulation precisely we need the following result (see [19] . Theorem 13.1). 

Theorem 1.8. For every <p G f) there exists (uniquely up to a complex phase) a coherent 
unitary : T — > T such that 

V$a(f)U+ = a(f) + (f,<l>)for allfet). 

Proof. We can proceed similarly the proof of Theorem 11.21 (see the Appendix) by translating 
the orthonormal basis {n^, n iM ) = (n ix \...ni M \) a*{ui M ) ni M ...a*(M il ) ni i |0) to 

U \n h , ...,n iM ) = (n il l..n iM \y l/2 [a*{u iM ) + {(f>,u iM )] ni M ...[a*(u h ) + {(f>, <OpiU |0) 
with the new vacuum 



|0) = exp 



exp[-a*(0)]|O> 



□ 



Remark. (i) The condensate vector <f> G fj needs not be normalized. Any pure state |^) (^>\ 
with ^ = |0) G T for some <fi G f) is called a coherent state. 

(ii) For generalized annihilation operators we get 

U;A(F)U^ = A(F) + (F, <P © J0)„er for all F G fj © If- 
Now we can describe the Bogoliubov variational states. Denote 

g" := {( 7 ,a,)|r w = ( l+ a Jir ) > 0,TV (7 ) < oo} . 

The Bogoliubov variational state p lja ,<f> associated with (7, a, <fi) G Q Bo x f) is defined by 

:= ^(IT^) for all 5 G 

where p 7)Q is the quasi-free state with the 1-pdm r 7jQ ,. In particular, the particle number 
expectation of the Bogoliubov variational state p lt a,<f> is 

p w (AO = Tr( 7 ) + ||0|| 2 . 

For a given Hamiltonian H : J 7 — > T and A > we can define the Bogoliubov ground 
state energy 

E B (\) = inf {p w (H)|( 7 ,a,0) G £ B x f),Tr( 7 ) + ||0|| 2 = A} 

where p lt(XjP (H) is the Bogoliubov energy functional and A stands for the total particle number 
of the system. 
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Remark. (i) Due to the variational principle, the Bogoliubov ground state energy E®(\) 
is always an upper bound to the quantum grand canonical energy 

E g (\) = inf{(tf,Mtf) : G F, ||*|| = 1, = A}. 

(ii) If iV G N then the grand canonical energy E S (N) is always a lower bound to the 
canonical energy 

N 

E(N) = inf{(*,ro) : * G (g), ||*|| = 1}. 

sym 

Moreover, if E(N) is convex w.r.t. N then E g (N) = E(N) for all N. 
Example 1.9 (A toy model). Let f) = M. and the Hamiltonian 

U = a*(l)a*(l)a(l)a(l). 
A straightforward computation shows that for A G N then the quantum energy is 

E g (N) = E(N) = N 2 — N 

and the Bogoliubov energy is 

E B (N) = inf \x 2 + x(4A - 2y/\(l + A)) + 2A 2 + A(l + A) 

\>0,x>0,\+x=N L 



= inf 

0<X<N 



N 2 + 2A(A - a/A(1 + A)) + A + 2Aa/A(1 + A) 
= N 2 -N + 0(N 2/3 ) as ->■ oo. 



Of our particular interest is the Bogoliubov variational theory for interacting Bose gases 
which we shall describe briefly below. 

Let f) = L 2 (Q) for some measure space Q with the inner product 



(u,v) = / u(x)v(x)dx. 
Jq 



In this case the mapping J : f) — > f)* is simply the complex conjugate, i.e. Ju(x) = u(x). 
Therefore, for simplicity we shall use notation 7 = J7J* and a = JaJ. 

The Hamiltonian consists of a one-body kinetic operator T, which is a self-adjoint operator 
on [), and a two-body potential operator W which is the multiplication operator corresponding 
to the funtion W(x,y) : f2 x — > R satisfying W{x, y) = W(y,x). The grand canonical 
Hamiltonian EI : T — > T can be represented in the second quantization as 



00 / N 
N=0 \i=l l<i<j<N 

= ^2(u m ,Tu n ) t) a* m a n + ^ {u m ® u n ,Wu p (g> u q )t, xt ,a* m a* n a p a q 



2 

m,n,p,q 



where a n := a(u n ) and {un}^^ is an orthonormal basis for () (the sum is independent of the 
choice of {u n }). 
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To represent the Bogoliubov energy functional explicitly in terms of (7, a, 0), it is conve- 
nient to introduce the integral kernel a(x,y) of the Hilbert-Schmidt operator which satisfies 



af)(x) = / a(x, y)f{y)dy for all / e L 2 (Q). 



n 

Similarly, we have the kernel 7(2;, y) of 7 and the density functional is formally defined by 
Py(x) := j(x,x). More precisely, because 7 is a positive semi-definite trace class operator, 
we have the spectral decomposition 7 = \ u i) ( u i\ anc ^ then we can define j(x,y) : = 

YjikUi{x)ui(y) and p^(x) := J2i ti\ui(x)\ 2 . Note that / p(x)dx = Tr(^y). 
Using the coherent transformations 

U^a n f/0 = a n + (u n , 0) 

and Wick's Theorem we find that 

Pw (H) = Tr(T7) + D(p^)+X(7,7)+X(a,a) 



+ Re 

hxn 



7(x, y)(j)(x)(j)(y) + a(x, y)(f>(x)</>(y) W(x, y)dxdy 



where 7 := 7 + \<f) (0|, Pj(x) = -y(x, x) and 



1 



D (f,g) = g // f( x )9(yW(x,y)dxdy, X(j, 7') = - / / 7(ar,y)7'(ar,j/)W r (x,j/)da:dj/. 



Here are some specific examples with respect to three cases: W > 0, W changes sign, 
and W < 0. 

Example 1.10 (Bosonic atoms). In this case we have 

f) = L 2 (M 3 ), T = -A-^-, W{x,y) 



\x\ \x — y\ 

We shall investigate the Bogoliubov theory for bosonic atoms in details in the next sections. 
In particular, we can show that the Bogoliubov ground state energy and the full quantum 
mechanics energy agree up to the leading order, and we conjecture that they even agree up 
to the second order. 

Example 1.11 (Two-component Bose gases). This is the case when 



f) = L 2 (M 3 x {±1}), T = -A x , W(x, e, y, e') 



ee 



\x-y\ 



It is already known that the Bogoliubov theory is also correct to the full quantum theory 
up to the leading order. More precisely, for large N, the correct leading term — AN 7 ^ 
was predicted by Dyson [7] using the Bogoliubov principle and then it was mathematically 
established by Lieb-Solovej [14] (lower bound) and Solovej [TS] (upper bound). 
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Example 1.12 (Bosonic stars). The system now corresponds to 



fj = L 2 (M 3 ), T = V-A + m 2 - m, W(x, y) = —. r 

\x-y\ 

where m > is the neutron mass and k = Gm 2 > 0. Up to the leading order, the ground 
state energy is approximated by the Hartree model [15]. Because the Hartree ground state 
energy is strictly concave, replacing the canonical setting by the grand canonical setting 
would make the energy much lower. Therefore, it is easy to see that the Bogoliubov ground 
state energy is much lower than the one of the full quantum model, although by adapting 
the ideas in jlOj we can show that the Bogoliubov variational model still has minimizers. 

2 Bosonic atoms 
2.1 Introduction 

For a bosonic atom we mean a system including a nucleus fixed at the origin in ]R 3 with 
nucleus charge Z > and N "bosonic electrons" with charge —1. The system is described 
by the Hamiltonian 

1=1 v 117 l<i<j<N 11 n 

acting on the symmetric space T-Ln = (8) S ym -^ 2 (^ 3 )- The ground state energy of the system 
is given by 

E(N,Z) =mf{(V,H NtZ V)\y en N ,\\*\\v> = !}• 

In fact, the ground state energy E(N, Z) does not change if we replace the symmetric sub- 
space n N by the full iV-particle space <g) N L 2 (R 3 ) = L 2 (R 3N ) (see, e.g., p3] p. 59-60). For 
usual atoms (with fermionic electrons), the Hamiltonian H^.z acts on the anti-symmetric 

N 

subspace /\(L 2 (IR 3 ) ® C 2 ) instead. For simplicity, we only consider the spinless electrons 

i=l 

because the spin number play no role in the mathematical analysis here. 

We recall some well-known fact about the full quantum problem. Due to the HVZ 
Theorem (see e.g. p3] Lemma 12.1), E(N, Z) < E(N - 1, Z) and if E(N, Z) < E(N - 1, Z) 
then E(N, Z) is an isolated eigenvalue of z- Unlike the asymptotic neutrality of fermionic 
atoms, in the bosonic case, the binding E(N, Z) < E(N — 1, Z) holds for all < iV < N C (Z) 
with ]im Z -+oo N C (Z)/Z = t c w 1.21 (see |Sl Oil El H] ) ■ 

The leading term of the ground state energy E(N, Z) is given by the Hartree theory [6]. 
In the Hartree theory, the ground state energy is 

E R (N, Z) = inf {8 u {u, Z) : | \u\ | 2 2 = N} 

where 

it . _ f .„ . Xl2j f Z\u(x)\ 2 f f \u(x)\ 2 \u(y)\ 2 , , 
£ H (u,Z)= / \Vu(x)\ 2 dx- / 1 . Y + // 1 V / M Y 1 dxdy. 

J J \ x \ J J \ x — y\ 

M 3 R 3 R 3 xR 3 
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By the scaling u{x) = Z 2 ui(Zx) we have 

S u ( Ul Z) = Z 3 S u ( Ul ,l). 

Therefore, 

E U (N, Z) = Z 3 e(N/Z, 1) where e(t) = E R (t, 1). 

It is well-known (see e.g. [5] [TT]) that e(t) is convex, e(t)' < when t < t c m 1.21 and 
e'(t) = when t > t c . Moreover, for any < t < t c ~ 1.21, e(t) has a unique minimizer <ft t , 
which is positive, radially-symmetric and it is the unique solution to the nonlinear equation 
h t <p t — where 

^ = -A-^ + K| 2 *^-e'(t). 
\x\ \x\ 

As a consequence, h t > 0. Moreover, since a css (h t ) = [— e'(t),0], there is a gap A 4 > if 
t < t c such that (/i t — A^P^ > where P^ t = 1 — P t with P 4 being the one-dimensional 
projection onto Span{0 t }. 

By scaling back, we conclude that E R (tZ, Z) has the unique minimizer and the operator 

ht,z = -A — -j— r + |0t,z| 2 * A - ZV W 
Fl Fl 

satisfies h tjZ 4>t,z = and ^ — Z 2 A t )P^ tz > when t < t c . 

Our aim is to investigate the first correction to the ground state energy E(tZ, Z). We shall 
analyze the Bogoliubov variational model for bosonic atoms and compare to the full quantum 
theory. From the general discussion on the Bogoliubov theory, we have the Bogoliubov 
variational problem 

E B (N,Z) = inf{£ B ( 7 ,a,0,Z)|( 7 ,a,0)G^ B xL 2 (R 3 ),Tr( 7 ) + ||0|| 2 = iV} (12) 

where 

£ B ( 7 ,a,0,Z) = Tr(-[A-Z|a;r 1 ]7) + J D(p^P^)+X(7,7)+X(a,a) 

+ [[ ^MM dxdy +Re ff ^ySEM dxdy . 

J J f — v\ J J f — y\ 

M 3 xR 3 M 3 xK 3 

Here we are using the notations 7 := 7 + \<f>) U>\ and 



2 J J \x — y\ 2 J J \x - y\ 

R 3 xR 3 R 3 xR 3 

The properties of the Bogoliubov theory for bosonic atoms are the following, which will 
be proved in the next subsections. 

Theorem 2.1 (Existence of minimizers). Let the nucleus charge Z and the electron number 
N be any positive numbers (not necessarily integers). 

(i) If the binding inequality 

E B (N, Z) < E B (N' , Z) for all < N 1 < N 

holds then E B (N,Z) has a minimizer. 
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(ii) The energy E B (N,Z) is strictly decreasing on N e [0, N C (Z)] with N C (Z) > Z for all 
Z and 

taint ^> (o 

Z^oo Z 



1.21. 



Theorem 2.2 (Bogoliubov ground state energy). If Z — > oo and AT/Z — t & (0,t c ) t/jen 

E B (N, Z) = Z 3 e(t) + Z 2 /i(t) + o(Z 2 ) 

[ 7 (x,y) + a(x,y)](f) t (x)(j) t (y) 



ait) := inf 

(7,")ee B 



Tr[ha] + Re 



f - y\ 



-dxdy 



The coefficient /i(t) is finite and satisfies the lower bound 

fi(t) < * -1 e(*) - e'(t) + /!(*) 

[l'(x, y) + y)}<f>t{x)<f> t {y) 



ait) := min 

( 7 ',a')GG s , 7 '0t=O 



Tr[/i t7 '] + Re J J 



\x-y\ 



< 0. 



Remark. (i) If we restrict the Hamiltonian H N Z into the class of TV-particle product func- 
tions ^ u = u <S> u <S> ... <8> m then by scaling = (AT — l)~ 1 / 2 Z 2 u (Zx) we have 



inf {V u ,H Ni zVu) = 

\u =1 



inf £ H (w , 1) = — re 



AT 



A" — 1 \\u\\ Q =(N ~1) / z N-l 
= Z 3 e(t) + Z\r x eit) - e'(t)} + o(Z 2 ). 

Because fj,(t) < t~ l e(i) — e'(t), the Bogoliubov ground state energy is strictly lower 
than the lowest energy of the product wave functions at the second oder. 

(ii) We believe, but do not have a rigorous proof, that the identity fi(t) = t~ 1 e(t) — e'(t) + 
Jl(t) holds and a minimizing sequence of /x(t) is given by 



7 = A 



110*11/ W 



+ 7, a = - V / A(l + A) 



<t>t 



110*11/ \\\<t> 



+ a' 



with A — > oo, where (7', a') is a minimizer for Jl(t). In fact, the upper bound 
li(t) > t~ l e{t) — e'(t)+]l(t) follows from the heuristic discussion on comparison between 
Bogoliubov energy and quantum energy below. 

We conjecture that the Bogoliubov theory determines the first correction to the quantum 
energy E(N, Z). 

Conjecture 2.3 (First correction to the leading energy). // Z — > 00 and N/Z = t e (0,t c ) 
then 

E(N, Z) = E B (N, Z) + o{Z 2 ) = Z 3 e(t) + Z 2 fi(t) + o{Z 2 ). 

A heuristic discussion supporting the conjecture is made in the last subsection of the 
article. While the picture is rather clear, some technical work is still needed to make the 
argument rigorous. 
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2.2 Existence of Bogoliubov minimizers 

To prove the first claim of Theorem 12.11 we shall follow the extending variational argument 
(see e.g. [12], Theorem 11.12). Before studying the variational problem E B (N,Z) in (fT2]h 
we start by considering the extended problem with the constraint Tr(7) < N, namely 

E B (<N, Z) = inf {£ B ( 7 , a, 0, Z) | ( 7 , a, 0) G G B , Tr( 7 ) + 1 |0| | 2 < iV} . (13) 

Lemma 2.4 (Extended problem). T/ie ground state energy E B (N, Z) is finite and decreas- 
ing on N. Moreover, the extended variational problem E B (<N,Z) in |73j) always has a 
minimizer. 

Proof. 1. By simply ignoring the non-negative two-body interaction and using the hydrogen 
bound, we have 



£ B ( 7 ,a,0,Z)>Tr 



A -) 

\x\ 



>-Tr(-A T )- — >-oo. (14) 



2. Next, we prove that E B (N\ Z) > E B (N, Z) for N' < N. For any trial state (7, a, 0) 
with (7, a) G £ B andTr7+||0|| 2 = N', choose 5 G C C °°(M 3 ) such that Tr( 7 ) + | |0| | 2 +| \g\ | 2 = JV 
and consider 

7 £ = 7 + |&) (&| 

where g e (x) = e 3 ^ 2 g(ex). Then (7 e , a) G £ B and Tr( 7e ) + ||0|| 2 = N. Moreover, since 
|V<7 e | — > in L 2 (R 3 ) and |g £ | 2 — >■ in L P (R 3 ) for any p > 1, a simple calculation shows that 

£ B ( 7s , a, 0, Z) £f( 7 , a, 0, Z) as £ -> 0. 

This ensures that E B (N', Z) > E B (N', Z). 

3. To show that E B (<N, Z) has a minimizer, let us take a minimizing sequence ( 7n , a n , n ) 
for E B (<N, Z). The lower bound (1T4"I) ensures that Tr(— A 7 „) is bounded. Consequently, 
all of ||7n(aJ,2/)||ffi(R3 xR a), ||an(z,y)||£rva(RS X R3) and ||0n||tfi( K 3) are bounded. By passing to 
a subsequence if necessary, we may assume that j n — ^ 7, a n — ^ 7, 0„ — ^ weakly in the 
corresponding Hilbert spaces, and their kernels converge pointwisely. It is straightforward 
to check that (7, a) G Q B and by Fatou's lemma, Tr(7) + ||0|| 2 < N. 

Fatou'lemma also implies that 



liminf Tr(— A 7 ) > Tr(— A 7 ). 

n— >oo 

The two-body interaction part of £ B (7„, a n , 0„, Z) can be rewritten as 

f f W( 7 n, «n, 0n) 

J J \x-y\ 

where 

Ty( 7n ,a n ,0 n ) = p ln (x)p Jn (y) + \>y n (x,y)\ 2 + \a n (x,y) + 0„(a;)0 n (?/)| 2 



+ 



P 7 n( a; )l0(z/)l + P-r n (y)\<P(?)\ +2Re(7 n (x,y)0 n (a;)0 n (2/)) 



> 0. 



Therefore, we may use Fatou'lemma again to obtain 

n^coJJ \x-y\ J J \x-y\ 
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Finally, because J~p^~ n — y/p^ in if^IR 3 ) we have the convergence 

Pin (g) , , /" PtM^ ^ 

— ; — : — ax — > / ———ax as n -)• oo. 



kr / x 



Therefore, we have 

liminf £ B (j n , a n , n , Z) > £ B (7, a, 0, Z) 



n— too 



and hence (7, a, 0) is a minimizer for E * (<N, Z). □ 

We now prove the existence of minimizers for the original problem E B (N, Z). 

Proof of Theorem\2J\ 1. If E B (N, Z) < E B (N', Z) for all < N' < N then any minimizer 
(7, a, 0) for the extended problem E B (<N, Z) must satisfy Tr(7) + ||0|| 2 = N, and hence it 
is a minimizer for E B (N, Z). 

2. That E(N,Z) is strictly decreasing on N £ [0, Z] follows by the same argument as 
in [5]. Assume that E B (N,Z) = E B (N',Z) for some < N' < N < Z. Let (7, a, 0) be a 
minimizer for E B (<N', Z). For any ip £ if 1 (M 3 ), let us consider the trial state (7 e , a, 0) with 

7 £ = 7 + e|^)(^|, £>0. 

For e > small we have Tr 7 £ + 1 10| | < iV and hence 

£ B ( 7e , a, 0, Z) > £ B (iV, Z) = E B (N', Z) = £ B ( 7 , a, 0, Z). 

Therefore, 

d 

+2ReX(7,|^)(^|). (15) 

On the other hand, let us replace cp by Pl(x) '■= L~ 3 ^ 2 p 1 (x/ L) where cpx £ if 1 (IR 3 ) such 
that pi is radially-symmetric and pi{x) = if \x\ < 1 and <^i(x) > if |x| > 1. Then for 
large L one has 

(<p L ,-A(p L ) = L~ 2 ((pt, —Aipi) = 0(L~ 2 ), 



tfe 



£ B ( 7£ ,a,0,Z) = ( v ,-Ap) L2 - [^^dx + 2D(p, / ,\p\ 2 ) 



x / \x 



Moreover, by Newton's theorem, 

max{|ar|,|2/|> 

R 3 xR 3 

and by Holder's inequality, 



dy, 



2 Re X(y, M M ) = //^l^« s 

J J \x-y\ 



>XJI 



/ \ 12 

< J J \l(x,y)\dxdy 
\\>L,\y\>L j 






Mv)\ 2 




x — 


y 


2 



1/2 

dxdy I = o(L~ l ) 



21 



Thus if we replace ip in ( TT5|) by tp L then we obtain 



< 0(U 



\Z - N')L 



A r-1 



which is a contradiction to the assumption N' < Z. Thus iV i— > E B (N, Z) is strictly decreas- 
ing when < N < Z. 

3. Now we show that E B (Z, N) is strictly decreasing on N G [Z, N C (Z)] with 

liminf N C (Z)/Z >t c ^ 1.21 

We shall need some properties of the Bogoliubov ground state in Lemma [221 which is derived 
in the next section. 

Take a large number Z and assume that N \— y E B (N,Z) is not strictly decreasing on 
Z < t'Z for a fixed value t' < t c . Then there exists N = tZ £ [Z, t'Z] and 5 > such that 
E B (N, Z) = E B (N + 5, Z) and E B (N, Z) has a ground state (7, a, 0). Because 



£ B ( 7 , a, 0, Z) = E B (N, Z) = E B {N + 5, Z) < £ B ( 7 , a, VTT^<p, Z) 
for e > small, we have 



< d 

~ de 



£ B (7,a,v / T + i0, Z) 



e=0 



-A - -r-r + ft * \.\ 1 



+ 



l{x,y)<f>(x)<f>(y) 



Re 



f - y\ 



a(x,y)<j)(x)<t>(y) 



\x - y\ 



Because 



we get 



Thus 



£ B ( 7 , a , 0) = E B (N, Z) < E H (N, Z) < Tr 



"A- FT J : 



n/(x,y)<j)(x)<f>(y) 



O<(0, -A- — 



\x - y\ 



1-1 



+ Re 



a{x,y)(f>{x)(f>{y) 



< 0. 



F - 1/1 



>A,z0> +2D(f*y- \4> t)Z \ 



e'(t)Z 2 ||0|f 



On the other hand, using the estimates in Lemma 12.71 we have 

{<f>,ht, z <f>) = o(Z 2 ), 
e\t)Z 2 \ |0| | 2 < e\t)Z 2 (tZ + o{Z)) = te\t)Z 3 + o(Z 3 ), 

D(p^-|0 iiZ | 2 ,|0| 2 ) < y/Dfa - \<j>t,z\ 2 ,Py - \^z\ 2 )-\Jd{\M 2 ^ I'M 2 ) = o(Z 5 / 2 ). 

Therefore, 

< (0, h t)Z 4>) + 2D(p^ - |0 tiZ | 2 , |0| 2 ) + e'(t)Z 2 \\<t>\\ 2 < te'(t)Z 3 + o(Z 3 ). 
However, it is a contradiction because te'(t) < when 1 < t < t' < t c . □ 
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2.3 Analysis of quadratic forms 



We consider the minimization problem /i(t) of the quadratic form in Theorem 12.21 Recall 
that 



//(*):= inf q t 0y, a) and /j,(t) := inf q t h',a') 



where 



Qt(l,ai) 



Tr[h tl ) + Re J J 



\x - y\ 



dxdy 



Lemma 2.5 (Analysis of the quadratic form gt(7, a)). For any < t < t c we have 

-oo < < r l e{t) - e'(t) + Jl(t). 
Moreover, the minimization problem Jl(t) has a minimizer (7', a') and Jl(t) < 0. 

Proof. 1. Because ^(7, a) is a quadratic form of (7, a), for considering the ground state 
energy we may restrict (7,0) into the class of quasi-free pure state, i.e. ota* = 7(1 + 7). 
Since 7 > is trace class and a T = a, we can write 

7(x,y) = ^ X n u n (x)u n (y), a(x,y) = - a/ A ra (l + A n )M n (x)M n (y), 

n n 

where A n > and {it n } n is an orthonormal family on L 2 (M 3 ). Then 

g t (7, a) = y^[A n (tt n , Mn) + A„] 

with 

= A n 



u n (x)u n (y)(j)t(x)(j) t (y) 



x - y\ 



y/X n (l + X n )Bje jj 



u n (x)u n (y)(f)t(x)4>t(y) 



x - y 



2. We may assume that X n (u n , h tt zu n ) + A n < for all n; otherwise, if X n (u n , ht,zu n ) + 
A n < then 

g t (7,a) > qtirf.of) 

where 

7' = 7 - X n \u n ) (u n \ , a' = a + a/A„(1 + A„) \u n ) (u^\ . 

We have the gap (u n ,h t u n ) > AtWP^UnW 2 for all n. Moreover, \Ke(D(u, u))\ < D(u,u) 
for all functions u we have 



(16) 



){X n - ^/X n {l + X n ))>2D (u n (j>t, Un<t>t) m &x < — -, — yX~n > . 
Thus it follows from the assumption X n (u n , h t z u n ) + A n < that 



A^A^IIP^II 4 < 4\D(u n (f)t,u n (j) t )\ 2 for all n. 



(17) 
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On the other hand, observe that 

\\P L u n \\ 2 + \\P L u m \\ 2 = 2- \\Pu n \\ 2 - ||Pw m || 2 > 1 for all m ^ n. 

Therefore, there exists (at most) an element %q such that \\P u n \\ 2 > 1/2 for all n ^ i . As 
a consequence, flTj) implies that 



J2 X n< 16A- 2 Munfau^l 2 < 4A- 2 J I 



Mx)\ 2 \My)\ 2 

\x - y\ 2 



dxdy < C. 



3. Using h t > and (I16p we have 
Qtil, a) > A io + y^ J Ai> -D{u io (f)t, u io <j) t ) - Y 2 VKD{u n (f) t , u n (p t ) 



> -D{u io 4)t,Ui^ 



1/2 



> 



1 



I0*(*)l 2 l&(v)l s 



1/2 



2 \ J J \x — y\ 2 

Pxl 3 



2A7 1 



A) I 2 

I0*(x)| 2 |0t(?/)| 2 



\x - 



> -C. 



4. To see the upper bound on fi(t) let us consider the trial state 



7 = A 



+ 7', a = -v^A(l + A) 



+ a' 



where (7', a') G £ B such that 7'^ = 0. One has 

Kt) < Qt(l, a) = 2 (A - V / A(l + A)) Diu^t, u^t) + q t {i , a'). 
Taking the infimum over all (7', a') and letting A — > 00 we obtain 

H(t) < -r^d^l 2 , \<p t \ 2 ) + Ji{t) = t~ x e{t) - e\t) + //(*). 

5. Now we consider The above argument shows that if {(7^, c>4)}r£=i is a mini- 

mizing sequence for Ji(t) then Tr(7^) is bounded. Therefore, it follows from the standard 
compactness argument that Ji(t) has a minimizer. To see that Jl(t) < 0, let us consider 



7' = A' \u) (u\ , a' = -y/X'(l + X') \u) (u\ 

where u is a normalized real- valued function in L 2 (IR 3 ) such that (u,(f>t) = 0. Because 
D(u(f> t ,u4>t) > we have 

Kt) < q t (j', a') = X'(u, h t u) + 2 (x' - y/X'(l + A')) t^D^t, u(j) t ) < 

for some A' > small enough. □ 
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Remark. The analysis here works out for a more general setting. For example, if h is a 
positive semi-definite operator on L 2 (Q) with inf cr ess (h) > and W is a positive semi-definite 
Hilbert-Schmidt operator on L 2 (Q) with a real-valued kernel W(x,y) then 



inf J Trf/ry] + Re / / {l{x, y) + a(x, y))W(x, y)dxdy | > — oo. 

(7,a)6S B 



By scaling 4>t,z(%) = Z 2 cf) t (Zx) 1 j(x) = Z 3r y'(Zx, Zy), a(x) = Z 3 a'(Zx, Zy) we have 



inf q t z (7 «) = inf Z 2 q t ('j', a') = Z 2 fi(t) 



where 



q t , z ( 1 ',a') = Tr[h t ^'} + J J 



y) + a'(x, y))(j)t,z(x)<frt,z(y) 



l 3 xM 3 



\x - y\ 



dxdy. 



To prove Theorem 12 .2[ we need to consider some perturbation form of qt t z- 

Lemma 2.6 (Analysis of pertubative quadratic forms). Let G L 2 (IR 3 ) such that \\<p\\ < 
||0t,z||; ||V0|| < CZ 3 / 2 and||P J -0|| < C where P 1 - = 1— P with P being the one- dimensional 
projection onto <pt,z- Then for Z large we have 



inf q t ,z{ , y,a,<f)) > 
(7,a)ee B 



\\M\ 



■z 2 fi(t) - cz 2 - 1 / 10 



where 



qt,z(l, a, 0) = Tr[/i tiZ 7] + J J 



l(x,y)(f>{x)(f>{y) f f a(x,y)(j)(x)(f)(y) 



\x - y\ 



\x - y\ 



Proof. 1. We first consider the case when Tr7 is small. Assume that Tr7 < Z 1 ! 2 £ , where 
e — 1/10. In the integral involved with 7, we use the decomposition 



<Kx)<f>(y) = P<P{x)P<P{y) + P0(x)P x 0(y) + P^(x)<f>(y). 
Observe that all terms involved with P^-cf) have negligible contribution. For example, 



<y(x,y)P ± (j)(x)(j)(y) 



R 3 xR 3 



\x - y\ 



< 2Tr(7 2 ) 1/2 ||P ± 0ll-l|V0|| < CZ 2 ~ £ . 



Thus 



le(x } y)(p(x)(p(y) 
\x - y\ 



> 



>y £ (x,y)P(j)(x)P(f)(y) 

\x - y\ 



-CZ 



2-e 



R 3 xR 3 



Together with the similar bound on the integral involved with a, we arrive at 
<M(7, a, 4) > (l- S Tr[hzl] + ir^QtAl, «) - CZ 2 ^ 



Ut,z\ 



> 



\\P<P\\ 2 

ll<MI 2 



Z 2 /j(t) - cz 2 - £ . 
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2. Now we consider the case when Tr7 is large. Assume Tr-y > Z 1 ! 2 £ . Following the 
proof of Lemma 12.51 we may assume that 



7 = Ai | ui) (ui\ + l', ol = -a/Ai(1 + Ai) \ui) (ui | + ol 
where = 1 and (7', a') is the 1-pdm of a pure quasi-free state such that 

TrY < C, AillP^MiH 2 < C and im = = aV 
Because Ai = Tr7 — Tr7 ; > Z 1 ^ 2 ^ 6 — C and AidP^iiiH 2 < C, we have 

\\p L Ul \\ 2 <cz- l i 2+£ . 

As a consequence, 

Tr(P 7 ') = ^A n ||P« n || 2 < Tr 7 '^||Pu n || 2 < Tr 7 '||PV|| 2 < Z~ 1/2+e . 

3. We shall compare qt,z(j, with %z{l" where 
-f = X 1 P\u 1 )(u 1 \P + P ± 'YP- L , 



a 



" = -y/X^l + XjP \ Ul ) (UT| P + PVP J 



It is easy to see that (7", a") G Q B . 

We first consider the terms involved with u\. We have 



Ai // - v/ATlTTAO Re 



k - y\ 



> (Ai - ^1(1 + AO) 

B 

Then we use the decomposition 



u 1 {x)u 1 {y)((){x)<i){y) 



F - 2/1 



ui(x)mi(?/) = Pux{x)Pu\{y) + Pui(x)P- L Mi(y) + P 
0M0(y) = P0(x)P0(y) + P<j>{x)P^<t>{y) + P^0(x)0(y). 

Note that all terms involved with either P u\ or P^-cf) have negligible contribution. For 
example, we have 



Pm 1 (x)Pm 1 (|/)P ± 0(x)0(|/) 



t 3 xt 3 



F - 2/1 



P ± m 1 (o;)m 1 (|/)P0(x)P0(|/) 



B 3 xR 3 



F - 2/1 



<2||P Ml || 2 ||P ± 0||.||VP0||<CZ 3 / 2 , 



< 2||P- L u 1 ||.|K||.||P0||.||VP0|| < CZ 2 - 1 / 4 ^/ 2 . 



26 



Thus 

ui(x)ui(j/)0(x)0(y) , . r p /" /" Mi(a:)ui(y)0(x)0(y) 



Al // - A(l + AQ Re || 



|ac - ?/| 



Mi(x)^i(y)0(x)0(y) 
|x - 2/| 



> (Ai - V M1 + AO ) / / PM*)Pui(vW(*)P<l>(v) _ CZ 2-i,w {u 



Next, consider the terms involved with (7', a'). In the integral, 



\x — y\ 

R 3 xR 3 

we use the decomposition 

y = P ± 1 'P ± + P^'P + P^', 

W)<P{y) = pW)P<P{y) + P~ I W)P^y) + P~ I W)<P{y)- 

Observe that all terms involved with either Py or P^ have negligible contribution. For 
example, we have 



(Py)(x,|/)P0(x)P0(2/) 



B 3 xR 3 



|x - y\ 



< 2[Tr(P(y) 2 P)] 1 /2||P0||.|| V P0|| 

< 2[Tr(y)] 1/2 [Tr(P 7 / )] 1/2 ||P0||.||VP(/.|| < CZ 2 - 1/8+e/A 



and 

(P ± 1 'P ± )(x,y)P I ^ix) ( f>(y) 



\x - y\ 



< 2[Tr(P ± (y) 2 P ± )] 1/2 ||P ± 0ll-l|VP0|| < CZ 3/2 . 



Thus 



y(x,2/)0(x)0(2/) > rr {p^jp^jx^p^pm _ cz ^ 1/8+£/ ^ (19) 

k - 2/1 ~JJ k - 2/1 



Similarly we have 



a'(x, y)<P(x)<P(y) > /" t (gV 'pJ-)(x,y)P<j>(x)P<i>(y) _ CZ 2-i/s+e/4_ (2Q) 



a? — 2/1 I ^ — 2/| 



Putting (TT8T) . (|T9p . (jzOjl together and using the fact h tt z > and ht,zP = 0, we obtain 

Il^ll IIWI 
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4. In summary, from Case 1 and Case 2 we have in any case 

Choosing e = 1/10 we obtain 

inf B g,z(7,«,0) > J^C^W " CZ^-Vio. 

(7,a)e5 B 110^11 



2.4 Bogoliubov ground state energy 



□ 



We are now ready to give the proof of Theorem 12.21 

Proof. Upper bound. Fix e > small. Choose (jt,e-> a t,e) £ such that 

Qt(%e, at ttE ) < pit) +s. 

Choosing7(x, y) = Z 3 ^ £ {Zx,Zy), a(x,y) = Z 3 a tt£ (Zx, Zy) and <p(x) = Z 2 <pt-Tr( lt , E )/z{Zx), 
we have Tr(7) + ||0|| 2 = tZ and 

£ B ( 7 , a, <P, Z) = Z 3 £f =1 (^_ Tr(7t , e)/z ) + Z 2 [q t ( lt>£ , a t>£ ) + e'(t) Tr( 7ti£ )] 
+Z [-D(p 7t , e ,p 7t , e ) +X(7 t)£ ,7 t ) + X(a t|E , a t , e )] 
< Z 3 e(t - Tr( 7t , £ )/Z) + Z 2 [/i(t) + e'(t) Tr( 7t , £ )] + ZC £ 
= Z 3 [e(t) - (Tr( Tt)£ )/Z)e'(t) + o^- 1 ) Tr 7e ] 

+Z 2 [//(*) + e + Tr( 7 ^)e'(f)] + ZC £ 
= Z 3 e(t) + Z 2 (p(t)+e + o(l)C £ ). 

Thus 

E B (N, Z) < Z 3 e(t) + Z 2 {p{t) +e + o(l)C £ )- 
Because e > can be chosen as small as we want, we can conclude that 

E B (N, Z) < Z 3 e(t) + Z 2 p(t) + o(Z 2 ). 

Lower bound. It suffices to consider (7, a, <ft) such that £ B (7,a,0, Z) < Z 3 e(t), and 
hence Tr[— A7] < CZ 3 . We shall denote by P the one-dimensional projection onto the 
Hartree ground state 4> tt z and P L = 1 — P. 

In the expression of £ B (7, a, 0, Z), if we ignore the non-negative terms X( 7 , 7), X(a, a) 
and estimate the direct term by 

D(pzj, ft) = 2D(p 7 , |0^| 2 ) - D(\^ z \ 2 , \^ z \ 2 ) + D(p^ - \<f) ttZ \ 2 , p 7 - \^ z \ 2 ) 

> 2D(p^, \^ z \ 2 ) - D(\<f> t;Z \ 2 , \4> t<z \ 2 ) = 2D(p^, \4> t<z \ 2 ) + Z 3 e(t) - Z 2 e\t) Tttf) 

then we arrive at 



£ B ( 7 , a, 0, Z) > Z 3 e(t) + Tr(^ ,7) + / / + Re / f ^ifl^M , (21) 



28 



By the same argument of the proof of Lemma 12.51 we have 

- Tr[ht zj]+ [f ^ y)W)m +Re[[ a ^^ x ]^ > - C Z 2 . 
2 ' JJ \x~y\ JJ \x-y\ 

R 3 xR 3 R 3 xR 3 

Putting this bound together with the gap Tr[h t 3-7] > A t Z 2 Tr(P- L 7) into ( 121]) . and comparing 
with the upper bound £ B (7,a,0, Z) < Z 3 e(t) we obtain | | J P- L 0[ | < C. 
We are now able to apply Lemma 12.61 to conclude from f[2"Tj) that 

P B ( 7 , a, 0, Z) > Z 3 e(t) + (0, h^) + j^Z^(t) - CZ 2 ^ . 

\Wt,z\\ 

Because ||-P</>|| 2 < ||</>t,z|| 2 = tZ, we obtain the desired lower bound. □ 

From the above proof of the lower bound, we also obtain the following estimates on the 
ground state, which will be useful in the proof of the binding up to the critical number t c Z. 

Lemma 2.7 (Properties of Bogoliubov minimizers). If (7, a, 0) is a minimizer for E B (N, Z) 
(or more generally, if P B (7,a,0,Z) = Z 3 e(t) + Z 2 p(t) + o(Z 2 )) then Trfp^) < C, 
(4>, h t ,z4>) = o{Z 2 ) and 

D(p^-\^ z \ 2 ,p^-\^ z \ 2 ) = o(Z 2 ). 

In particular, it follows from (<j),h t .z<p) = o(Z 2 ) that ||-P0|| 2 = tZ + o(Z). Here P is the 
one- dimensional projection onto the Hartree ground state <pt,z- 



2.5 Comparison to quantum energy: a heuristic discussion 

Let us discuss on the comparison between the Bogoliubov ground state energy E B (N,Z) 
and the quantum energy E(N, Z) in Conjecture 12.31 

First at all, due to the variational principle, the Bogoliubov energy E B (N, Z) is a rigorous 
upper bound to the quantum grand canonical energy 

00 

E%N, Z) = inf{(*, H NjZ ^f), # e T, 11*11 = 1}. 
jv=o 

It is believed that the ground state energy E(N,Z) is a convex function on N (see |13j . 
p. 229), which is equivalent to E g (N,Z) = E(N,Z). If this conjecture is correct then the 
Bogoliubov energy E B (N, Z) is also an upper bound to the canonical energy E(N, Z). 

In the following, we shall argue heuristically why the Bogoliubov energy E B (N,Z) is a 
lower bound to E(N, Z) (up to an error o(Z 2 )). Some further work is required to make the 
argument rigorous. 

Choosing an orthonormal basis {u n }^ =0 for f) with u = 4>t,z/\\<f>t,z\\, we can represent 
the Hamiltonian Mz = ©jv=o H^,z in the second quantization 

Hz = h myn a* m a n + - W TOi „ )P)? a^a*a p a 9 

m,n>0 m,n,p,q>0 
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where a n = a{u n ) and 



hm,n ("mi ( A Z\x\ )lt n ), ^Vm,n,p,q 



u m (x)u n (y)u p (x)u q (y) 



\x - y\ 



Assume that \I> is a ground state for E(N,Z). We shall denote by (IHhA^ the expectation 

Step 1. As in [3] we have the condensation Tr(P ± 7^ I ) < C where P is the one- 
dimensional projection onto uq. Let us denote 7 = P ± 7*P- L , a = P^o^P -1 , Ao = a^ao 
and A = (JVoV Then (7, a) e £ B and A - A = Tr( 7 ) < C. 

Step 2. The leading term Z 3 e(t) of the ground state energy E(N, Z) comes from the 
terms of full condensation, namely /ioo a o a o an d W / oooo a o a o a o a o- Similarly to the computation 
to the energy of product functions, we have 

h 00 (cto a o)* + W^oooo (ao a o a o a o)^ 
= (uq, (-A - Z\x\- 1 ) u ) N + ((A^ - No)D(\u \ 2 , KH 

> (uq, (-A - Z\x\- 1 ) u ) N + (N^ - No)D(\u \ 2 , \u \ 2 ) 

^ N Z 3 {No -I 

> re 



JV -1 \ Z 

= Z 3 e(t)-Z 2 e'(t)Tr( 1 ) + Z 2 [r 1 e(t)-e'(t)} + o(Z 2 ). (22) 

As a consequence, the expectation of the rest of the Hamiltonian should be of order 
0(Z 2 ). 

Step 3. Because almost of particles live in the condensation m , we may hope to eliminate 
all terms W mnp the two-body interaction which have only or 1 operator a% 

(where af is either a or Oq). 

Step 4. Now we apply the Bogoliubov principle in which we replace any af by \/No ~ 
y/N. We can see that the terms with 1 and 3 operators af should be canceled together. In 
fact, 

{hom{^o a m) ^ + WoOOm { a a a a m) 

m>l 

~ {hom(a>O a m) $ + W / OOOm(OoO m )^) 

m>l 

= ^ ("m, (-A - + A|w | * M ) (ao a m)^ = 

m>l 
m>l 

Here we use the fact that Uq is the ground state for the Hartree mean-field operator 

ht tZ = -A - Z\x\- 1 + |<M 2 * l.r 1 - Z 2 e\t) 
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It remains the terms with precisely or 2 operators af, 

(h mn (a* m a n ) y + W m QQ n (a* m a* aQa n ) ^) 

m,n>l 
m,n>l 

= Tr [(-A - Z\x\- 1 + iV|u | 2 * I.]" 1 ) 7] (23) 

and 

) 9 + Re[W mn00 ( 

a m a n a O a o) *]) 

m,n>l 

{ a m a n)^j + ^R- e [W / mnOo( a m a n)*]) 

m,n>l 

= Re // + °(*>v))M*)<i>tAv) dxd y. (24) 

R 3 xR 3 

Step 5. Putting the approximations (1221 . (123"|) and (12I]) together we obtain the desired 
lower bound 

(Hz)* > ZV^ + ^Ir^-e'^] 

+ Tr[/^ 7 ] + Re // ^ + ^>^W^) + o(Z ^ 

JJ \ x ~y\ 

R 3 xR 3 

Because (7, a) G {? B and 70t,z = one has 

^[r 1 ^) - e'(*)] + Tr[/^ 7 ] + Re// + ^g^M^ > Z ^(t). 

R 3 xR 3 

Thus we arrive at the desired lower bound 

'3 „//_.\ 1 72,,/A 1 / ry2\ 



h)*>Z 6 e'(t) + Z 2 fi(t) + o(Z 2 ). 



Appendix 



Proof of Lemma ll.li It is obvious that T > if and only if 7 > 0, a* = JaJ and 
(/ © Jg, Tf © Jg) = (f, 7/) + (g, (1 + i)g) + 2 Re(a Jf, g) > 0, V/, g G I). 
Using a simple scaling g — tg, t G C, we can see that the latter inequality is equivalent to 

{f,lf){g,{l + l)g)>\{aJf,g)\\ V/^ef). 
Replacing g by (1 + , y)~ 1 g, we can rewrite the latter inequality as 

C/W/Xs, (1 + l)- l g) > \(aJf, (1 + 7)-^) | 2 , Vf,g G f). (25) 
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Note that (j2J) follows from (J25]) by choosing g = aJf. Reversely, we can see that (J2J) 
implies ( 12 5 p by using the Cauchy-Schwarz inequality for the positive definite quadratic form 
Q(u,v) = (u, (1 + 7)~ 1 w), i.e. 

(/, 1 f)(g, (1 + 7)^) > {aJf, (1 + X)~ 1 aJf)(g, (1 + 7)^) > | (a J/, (1 + 7)^)T • 

□ 

Proof of Theorem li.H The proof below follows [Tj5] Theorem 9.5 (sufficiency) and [TH] The- 
orem 6.1 (necessity). 

Sufficiency. Assume that VV* is trace class on f). We shall construct the unitary Uy- 

1. Let {ui}i>i be an orthonormal basis for fj. Recall that an orthonormal basis for 
^^'^(f)) is given by 

\n h ,...,n iM ) = (n il \...n iM \y 1/2 a*(u iM ) ni M ... a *(u h ) n ^ |0) , 

where rij run over 0, 1, 2, ... such that there are only finite rij > 0. 

We start by constructing the new vacuum |0) v = Uy |0) which is characterized by 

A(V(w 4 ffiO))|0) v = 0. 

for all i = 1,2,..., namely 

A(V{u t © 0)) = A{Uui © JVJui) = a(Uui) + a*(VJui) 

are the new annihilation operators. 

2. The first step is to choose an convenient basis {ui}. From V*SV = S = VSV* we have 

UU* = 1 + VV*, U*U = 1 + J*V*VJ 

and C = C* where C = U*VJ. Since U*U — 1 is trace class, U*U has an orthonormal 
eigenbasis on I). On the other hand, because U*U commutes with the conjugate linear map 
C = U*J*V and C*C is trace class on f}, we can find an orthonormal basis {«i}i>i for f) 
consisting of eigenvectors of U*U such that they are also eigenvectors of C. 

Denote fii := Ht/WiH > 1 and fi := fi~ 1 Uui. Then {fi}i>i is an orthonormal basis for f). 
Since 

(fj, VJui) = ^(Uuj, VJui) = ^~ x (uj, Cui) = 

for all j 7^ i, we must have VJu^ G Span{/j}. Note that if we change Mj's by complex phases 
then it still holds that (uj, Cui) = for all i 7^ j (although Wj's maybe no longer eigenvectors 
of C). Therefore, we can change Uj's by complex phases to obtain VJui = V{f for some 
Ui > 0. Thus there is an orthonormal basis {fi}i>i for f) such that the new annihilation 
operators are 

A(V( Ui © 0)) = tna(fi) + m*(fi), i = l,2,... 
where m > 1, v t > 0, fif - vf = 1 and J2i=i v l = ^(VV*) < 00. 
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3. This representation allows us to construct the new vacuum |0) v explicitly 



10), 



M oo , 

j5s.n(i-c^) a ) V4 E(-^ 

7=1 n=0 v ^ J 



,2n 



nl 



10} 



n(l-(^M) 2 ) V4 exp 

3=1 



10) 



It is straightforward to check that |0) v is well defined and is annihilated by the new annihi- 
lation operators A(V{ui © 0)). Having the new vacuum |0) v , we can define |n i:l , ...,n iM ) v = 
i]\n il ,...,n iM ) by 

\n h ,...,n iM ) v = (n h L.n iM \)- 1 / 2 A*(V(u iM © 0)) n ^ ...A*(V(u h © 0)) n n |0) v . 

4. Finally we need to prove that the new vectors |n ii; ...,n iM ) v indeed form a basis for 
J 7 . The trick is to use the formula 



io)=n( i -(^) 2 )" i/4ex p 



10), 



and express the old basis vectors [n^, ...,rii M ) in terms of the new ones. Since the new 
vectors {n^, n iM ) v span all of the old basis vectors {n^, ...,n iM ), the new ones span the 
whole space J 7 . 

Necessity. Assume that there exists a normalized vector |0) v G J 7 such that A(V(u © 
0)) |0) v = for all u 6 f). We shall prove that V^V™ must be trace class on f). 

5. Let |0) v = ®n=o^n where ^ e (g) s ^ m f). Then the condition A(V(« © 0)) |0) v = 
is equivalent to 



a(Uu)Vi = and a(Uu)V N+2 + a*(VJu)V N = for all u e f), JV = 0, 1, 2, .. 



(26) 



Since £/£/* = 1 + K^* > 1 we have Ker(C/*) = {0}, and hence Ran(C7) = fj. Therefore, 
it follows from a(Uu)^i = for all u 6 f) that = 0. Then, by induction using (|2"E)) we 
obtain *i = ^ 3 = vj/ 5 = ... = 0. 

If = then the same argument deduces \l/o = ^2 = ^4 = ■■• = which contradicts 
with |0) v ^ 0. Thus ^ e C\{0} and from ([26]) with iV = we have 

a(Uu)^ 2 + ^ VJu = for all u E f). (27) 

6. Introducing the conjugate linear map if : f) — )■ f) defined by 

v? 2 ) = (*2, <Pi ® ¥>a) for all v?i, ^2 e f). 

A straightforward computation shows that Ty(H*H) = || l I / 2 || 2 . Moreover using fl27|) and the 
symmetry of ^2 we have 

(-* VJ^i,^ 2 ) = (a(C/^i)# 2 ,v? 2 ) = (* 2 ,a*(^i)V2,) 
= V2(V2,U(p 1 ®(p 2 ) = V2(HU<p 1 ,<p 2 ) 

for all </?i, (^2 G f). This means — ^0^^ = Because Z7 is bounded and is trace class 
on f), we conclude that 



VV* = 2% 2 HUU*H* 



is trace class on f). 



□ 
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The rest part of proof of Theorem \l.b\ We now prove that V is the 1-pdm of the state p = 
Tr[G] _1 G. Recall that J 7 has the orthonormal basis 

Kn 2 ,...> = (n 1 W...)- l/2 (a* i r(a* 2 r...\0) 

where |0) is the vacuum and n\, n 2 ... run over 0, 1, 2, ... such that there are only finite rij > 0. 
A straightforward computation shows that 

Tr(G) = 22 ( n i> n 2, G|rii,n 2 , ...} 

n 3 =0,l,2,... 

(niW....)- 1 (0| II exphA^KD |0) 
rij=o,i,...;jei iei 

= £ (niM...)" 1 (01 II (V £ ( ~ A ° fc fc f Q ° fc (aI)-A |0) 
^=0,1, ■■■;ie/ te/ V fc=o ' ' / 

= e (-M...)-<oinff: ( " e ' ) V ( "'" > )i o > 

fii=0,l,...;j6l ig-T \fc=0 / 



e n^"=ni^ 



nj=0,l,...;j67 i iei 



t- < OO 



since J2iei e ^ < 00 • Thus p is well-defined. 

We check that T is indeed the 1-pdm of p. Note that \ni,ri2, ■■■) and G |ni,ri2, •••) have 
the same number of particle Ui for any i = 1, 2, .... By the same way of determining Tr(G) 
we find that Tr(ajOjG) = and 

Tr(a*a>jG) = 5^ Tr(a*ajG) 
= II + ( £(^ ! ) _1 (0|ar i a> l exp(-c l a> J )(a*rM0) 



oc 



*«( II (! + ^) ) ( E ("iO" 1 (01 £ V (a?)"* |0) 

%( n d+^l fEM-<Dif; ^y^ io) 

\kel,k& J \n z =0 r=0 

I U (1 + Afc) I I ^ exp(-Qni)ni J = <%Aj JJ (1 + A fc 
\fcei,fc^t / \ni=o / feel 



in which we have used 



oo , oo tZ 1 

^Texpt-c^K = ^ exp (- CiWi ) = -- 1 _ 

ri;=0 n.;=0 



eXp( - Ci) -A,(l + A, 



(1 - exp(-Cj)) : 
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From the above computations we find that 

p(ajOj) = (Tr(G)) -1 Tr(ajajG) = = {u^aJuj 

and 



p(a*aj) = (Tr(G)) 1 Tr(a*a,jG) = SijXi 



for any Thus T is indeed the 1-pdm of p. 

3. Finally, we check that p is a quasi-free state. One way to do it is to consider p as a 
limit of appropriate Gibbs states, see jl] (eq. (2b. 34)). In the following, we shall give a more 
direct approach by mimicking the proof of Wick's Theorem in [H] . 

It suffices to prove ©-(IZD when A(Fj) is either a creation or annihilation operator, which 
we denote by q. Our aim is to show that 

Tr[ciC 2 c 3 c 4 ...c fc G] = ~T^^p Tr [ c 3 c A ...c k G} (28) 

Tr[cic 3 G] Tr[cic fc G] 
Tr[G] Tr t C2C4 "- CfcG ] + -■ + Tr[G] T W3-- c fc-i G ] 

and the result follows immediately by a simple induction. By the same way of computating 
Tr[G] we may check that 

Tr[cic 2 G] 

Tr j G j = f(ci)[ci,c 2 ] (29) 
where [ci,C2] = cic 2 — c%c\ G {0, —1, 1} and 



[1-e A 1 if d = aj, j G I, 
[l-e^y 1 if a = 
L if d = dj, j i I, 
if ci = a*, j £ I. 



Thus (BHD i s equivalent to 



Tr[cic 2 c 3 C4...c fc G'] = /(ci)[ci, c 2 ]Tr[c 3 c 4 ...c fc G'] (31) 
+/(ci)[ci,c 3 ]Tr[c 2 c 4 ...c fc G'] + ... + f(c 1 )[c 1 ,c k \Tr[c 2 c 3 ...c k _ 1 G]. 

We can prove (131 j) as follows. From the identity 

CiC 2 C 3 C4...Cfc = [Ci, C 2 ]c 3 C4...Cfc + ... + C 2 C4...Cfc_i [Ci, Cfe] + C 2 C 3 C4...CfcCi 

we deduce that 

Tr [ciC 2 c 3 C4...c fe G] = Tr [[c u c 2 ]c 3 c 4 ...c k G] (32) 

+••• + Tr [c 2 C4...c fc _i[ci,c fc ]G ! ] + Tr [c 2 c 3 c 4 ...c fc ciG'] . 

We first consider when c\ is either Oj or a* with j G /. In this case it is straightforward 
to see that C\G = e j CxG where (+) if C\ = a* and (-) if c\ = aj. This implies that 

Tr [c 2 c 3 c i ...c k c 1 G] = e ±A ^Tr [c 2 c 3 c 4 ...c fc Gci] = e ±Aj Tr [c l c 2 c 3 c 4 ...c k G] . (33) 
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Substituting ( 13"3"j) into ( 13"3"j) we conclude that 



Tr [cic 2 c 3 C4...c fc G] = - — '-^-Tr [c 3 c 4 ...c k G] 



[ci,c 2 ] 
1 - e ±x > 

[ci,c 3 ] r [ci,c fc ] 

+ 3 zr Tr c 2 c 4 ...c fc G + ... + T3-Tr c 2 c 4 ...c fc _iG 

1 — e ±A J 1 — e ±A J 



which is precisely the desired identity f[3"Tj) . 
If ci = a, for some j ^ I then 



Tr[c 2 c 3 c 4 ...c fc CiG'] = 



since djG — and (|3"T|) follows from f[3"3"j) . 
Finally if cj = a* for some j £ I then 

Tr[cic 2 c 3 C4...c fc G] = Tr[c 2 c 3 c 4 ...c fc Gci] = 

since Ga* = and we obtain (l3Tj) . □ 

Acknowledgments: I thank my advisor Jan Philip Solovej for giving me the problem and 
various helpful discussions. 
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